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The model of two electrons with Coulomb interaction on a two-dimensional (2D) disordered 
lattice is considered. It is shown that the interaction can give a sharp transition to delocalized 
states in a way similar to the Anderson transition in 3D. The localized phase appears when 
the ratio of the Coulomb energy to the Fermi energy becomes larger than some critical value 
dependent on the disorder. The relation to the experiments on metal-insulator transition in 
2D is also discussed. 



1 Introduction 

According to D. TsuiS, "the important thing is the interplay between disorder and electron- 
electron interactions. The FQHE (fractional quantum Hall effect) is, in some sense, the clean 
limit. But there's another limit, where both interaction and disorder are important ... there's 
always some disorder." Indeed, the recent experimental discovery of metal-insulation transition 
in two dimensions (2D) by Kravchenko et al. i attracted a great interest to this problem. 
This transition is especially surprising since according to the well established theoretical result 
H all states of non-interacting electrons in 2D disordered potential are exponentially localized. 
However, in reality the electron-electron interaction is present and the original xesuItE ,_as well 
as the new results of different groups in experiments with different materials □Ij'BlZI'B'Emy, show 
that the interaction can induce metallic behavior. Indeed, the majority of experiments are done 
in the situation where the parameter = l/y^vfnja^ ~ E(.e/Ep ^ 1. Here, E^e is the energy 
of Coulomb interaction, Ep is the Fermi energy determined by the charge density Ug and a*^ 
is the effective Bohr radius. In some experiments the value was as large as 10 - 30. In this 
situation the electrons are located far from each other and in a first approximation it is natural 
to consider the problem of only two electrons with Coulomb interaction. The first consideration 
of two particles with strong attraction was done by Dorokhov but it was ignored by the 
community. The studies of two interacting particles with short range interaction showed that 



repulsive/attractive interaction can lead to a stroneincrease of localization length or even to 
delocalize pairs of particles in dimension d > According tollllll in 2D the pairs 

of particles remain localized and their localization length Ic grows smoothly with the increase 
of disorder strength U or one-particle localization length li: ln(/c/^i) ~ k > with k ~ ^2P2, 
where r2 ~ U'^/Vli is the interaction induced transition rate, p2 ~ ^i/^ is two-particle density 
of states in the middle of the band and V is the hopping strength proportional to the energy 
band size B (B = AV for weak disorder). The case of the long range Coulomb interaction 
requires separate analysis. Generally, one can expect that the delocalization effect will be even 
stronger in this case since the particles are always interacting, in a difference from a short range 
interaction case. 



2 Analytical estimates 

The first estimates for two electrons with Coulomb interaction on a 2D disordered Anderson 
lattice were presented in 113. The lattice is characterized by the nearby hopping V and the 
diagonal disorder in the interval [— VF/2, M^/2], while the interaction is U/\ri — r2|. Then the 
parameter = U / {2V yjims) where Ug is the filling factor. If the distance R between the 
electrons is much larger than the one-particle localization length li (In/i ~ {y/W)"^) then the 
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Figure 1: Dependence of the inverse participation ratio ^2 on rescaled one-particle energy e/B for 
U/V =2, L = 16 (r, = 6.38, B = AV) and W/V = 5, 7, 10, 15 (from up to down). 

two-body coupling appears ^ly in the dipole-dipole interaction term. This gives the typical 
matrix element Ug ~ U / and the transition rate r2 ~ UgP2, where still p2 is determined 
by the estimate given above since the electrons can have a jump only on a distance from the 
initial position (otherwise wave function overlap drops exponentially). As the result the two 

electron levels become mixed by interaction when /tg = Xe ~ ^2P2 ~ (^l^^/'^s)^ > where 
is the value of at the density Ug = l/lf (one electron in a box of li size). For Xe > 1 the 
Coulomb interaction leads to a delocalization of two electrons in a way similar to 3D Anderson 
transitior0. Indeed, in this case the hopping goes effectively in 3D: the center of mass moves in 
2D and in addition the electrons slowly rotate around it that gives 3 dimensions. The rotation 
goes on a ring of width li and of radius i? ~ ^ li (for [/ ~ F); the size of the ring is fixed 
by the energy conservation e ^ U/R. Due to that the length L changes sharply from Ic ~ li to 
Ic ~ /i exp(7r?]^^'^Ke) when Kg crosses the critical value Ke ~ It is interesting to note that, 
as in the experiments (see Refs. 2,4-10), the localized phase corresponds to the large values of 
Ts'. physically the two-body interaction becomes weaker at low density. The diffusion rate in 
the metallic phase can be estimated as Dg ~ ^1^2 ~ Vne/li- Near the critical point Kg ~ 1 the 
diffusion rate (conductivity) drops with the decrease of disorder (increase of ^i). These estimates 
are done for the excited states in the middle of the band. 

3 Numerical results 

The above problem of two electrons in the 2D Anderson model in the triplet state is studied 
numerically. The maximal lattice size is L = 24. The numerical diagonalization is done in the 
following way: the Hamiltonian is rewritten in the basis of noninteracting eigenstates, from which 
only first M low energy one-particle states (orbitals) are selected and after that the Hamiltonian 
is diagonalized exactly. The special check is done to ensure that the low energy states are not 
effected by the above cutoff (e.g. by changing M in few times). Usually ND=4000 disorder 
realizations are used to average the fluctuations. The fact that the effect of interaction strongly 
depends on li (or W) is demonstrated in Fig. 1. Indeed, here the number of noninteracting 
eigenstates ^2 contributing in an eigenfunction at fixed interaction U/V = 2 is increased in 
about 50 times only by the change of the disorder W . This confirms the analytical result 





Figure 2: Dependence of the rescaled transition rate r2 (defined via the relation r2 = (^2 — l)/p2(e) with 
p2(e) being the two-electron density of states) on e/B for the parameters of Fig.l and the same order of curves, 



according to which the effect of interaction becomes stronger for larger li since ^2 ~ ^2P2{^)- 
The last relation allows to determine numerically the dependence of the transition rate T2 on the 
excitation energy E = 2e, counted from the ground state. This dependence is presented in Fig. 
2 and for a moderate disorder shows a very flat dependence on e and even a certain increase of r2 
very close to the ground state. If to assume that the matrix element Ug is independent of e then 
r2 should drop linearly with e since P2(e) ~ Zfe/y^lll. However, for a short range interaction 
is has been shown that for localized states r2 can be independent of e due to enhanced return 
probability near the Fermi level 0. For the case of long range interaction similar effects can be 
responsible for the flat variation of r2 with e in Fig. 2. More detailed studies are required to 
understand the properties of r2 near the ground state. 

Another part of numerical studies is devoted to the investigation of the level spacing statistics 
in the above model. Indeed, it is known that the localized phase is characterized by the Poisson 
distribution Pp(s), the metallic phase has the Wigner-Dyson statistics Piyn while the critical 
transition point has an intermediate statistics independent of the system sizeEj. It is convenient 
to study the transition between two limits with the help of the parameter rj = Jq°{P{s) — 
PwDis))ds/ fQ°iPp{s) — PwDis))ds, where sq = 0.4729... is the intersection point of Pp{s) and 
Pwois)- In this way r] = 1 corresponds to Pp{s), and ry=0 to Pwois)- The dependence of 
r/ on e is determined in the following way. For each disorder realization the spacing between 
nearby energy levels Ei is determined and then is averaged over ND disorder realizations for 
each i giving the P{s) statistics and 77 as a function of averaged excitation energy e = E/2. At 
higher energies the values of rj are in addition averaged in a fixed energy interval. In this way 
the total statistics obtained for P{s) and r] varies from NS = 12000 for low energy states up to 
NS = 10^ at high energy with high density of levels. Un example of rj variation with energy 
for different system size L is shown in Fig. 3 (see also Fig. 1 inlll for a stronger disorder). 
For large L the statistics becomes close to Pp{s) at low energy and to Pwnis) at e larger than 
a critical energy €c dependent on the disorder and independent of L. The transition in the 
spectral statistics can be qualitatively understand on the basis of the estimates given in the 
previous section. Indeed, since the interaction energy is C//i? ~ e the high energy states allow 
to have particles closer to each other {R ~ U/ e) that increases their interaction and finally leads 
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Figure 3: Dependence of ?? on e/B for U/V = 2, W/V = 5 and different L: 10 (o); 12 (□); 16 (full diamond); 20 

(*); 24 (x); so that 2.39 < < 9.57. 

to delocalization for R > l'^^^{U ~ y). In agreement with this picture the critical energy Cc 
decreases with the increase of li (decrease of W/V) as it can be seen from Fig. 1 in EH and 
Fig. 3. The fact that an interaction increases the locahzation length for two particles in 2D 
has been also seen in the other numerical simulationsE^il. However, the claim made there that 
the short range interaction gives a transition from localized to delocalized states is in a sharp 
contradictionwith the theoretical argumentsllUllH and probably shouldbe attributed to small 
sizes used inE3'E2l. The numerical data for the Coulomb case presented inEl are somewhat similar 
to the data preservted here and inlll, even if any theoretical arguments in the favor of transition 
were presented inEil. 

The variation of P{s) with the interaction strength is shown in Fig. 4. At small U the 
statistics approaches to the Poisson distribution while with the increase of U it tends to the 
Wigner-Dyson case. In the vicinity of the critical point €c the statistics is close to the critical 
statistics in the 3D Anderson model with periodic boundary conditions (see Fig. 4). This 
gives one more support for the physical picture developed in the previous section according to 
which two electrons in 2D are delocalized in a way similar to the 3D Anderson transition. 

4 Conclusion 

The analytical and numerical results obtained show that the Coulomb interaction leads to 
delocalization of two electron states in a way similar to the Anderson transition in 3D for 
fL<fs< "i^t!^- The model is restricted only by two interacting electrons and has delocalization 
only for excited states that represents its weak point. However, it gives a picture qualitatively 
similar to experimentally observed metal-insulator transition in 2D i and therefore it can be 
useful for a future complete theory. 
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The full line shows the critical P{s) in 3D Anderson model {W/V — 16.5, L = 14, taken from Ref. 23); the dashed 

lines give Poisson statistics and Wigner surmise. 



3. E. Abrahams, P. W. Anderson, D. C. Licciardello, and T. V. Ramakrishnan, Phys. Rev. 
Lett. 42, 673 (1979). 

4. S. V. Kravchenko, D. Simonian, M. P. Sarachik, W. Mason, and J. E. Furneaux, Phys. 
Rev. Lett. 77, 4938 (1996). 

5. D.Popovic, A. B. Fowler, and S. Washburn, Phys. Rev. Lett. 79, 1543 (1997). 

6. P. T. Coleridge, R. L. Wilhams, Y. Feng, and P. Zawadzki, Phys. Rev. B 56, R12764 
(1997). 

7. Y. Hanein, D. Shahar, J. Yoon, C. C. Li, D. C. Tsui, and H. Shtrikman, Phys. Rev. B 
58, R7520 (1998). 

8. M. Y. Simmons, A. R. Hamilton, M. Pepper, E. H. Linfield, P. D. Rose, and D. A. Ritchie, 
Phys. Rev. Lett. 80, 1292 (1998). 

9. V. M. Pudalov, G. Brunthaler, A.Prinz, and G. Bauer, JETP Lett. 68, 534 (1998). 

10. V. Senz, U. Doetsch, U. Gennser, T. Ihn, T. Heinzel, K. Ensslin, R. Hartmann, D. Gruet- 
zmacher, cond-mat /9903367 . 

11. O. N. Dorokhov, Sov. Phys. JETP 71(2), 360 (1990). 

12. D. L. Shepelyansky, Phys. Rev. Lett. 73, 2607 (lQa4); p. 201 in@. 

13. Y. Imry, Europhys. Lett. 30, 405 (1995); p. 211 inEl. 

14. D. Weinmann, A. Miiller-Groeling, J.-L. Pichard, and K. Frahm, Phys. Rev. Lett. 75, 
1598 (1995); p.221 in0. 

15. F. von Oppen, T. Wettig, and J. Miiller, Phys. Rev. Lett. 76, 491 (1996); p. 235 inEl. 

16. Correlated fermions and transport in mesoscopic systems, Eds. T.Martin, G.Montambaux 
and J.Tran Thanh Van, Editions Frontieres, Gif-sur-Yvette, (1996). 

17. D. L. Shepelyansky, |cond-mat/9902246t 

18. Ph.Jacquod and D.L. Shepelyansky, Phys. Rev. Lett. 78, 4986 (1997). 

19. B. I. Shklovskii, B. Shapiro, B. R. Sears, P. Lambrianides and H. B. Shore, Phys. Rev. B 
47, 11487 (1993). 

20. M.Ortuno and E.Cuevas, Europhys. Lett. 46, 224 (1999). 

21. E.Cuevas, |cond-mat/980813g| v2. 

22. I. Kh. Zharekeshev, and B. Kramer, Phys. Rev. Lett. 77, 717 (1997). 

23. D. Braun, G.Montambaux and M.Pascaud, Phys. Rev. Lett. 81, 1062 (1998). 



